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Abstract 

We present results for the spectrum of bb bound states in the quenched 
approximation for three different values of the lattice spacing, in the range 
0.05fm to 0.15fm. We find our results for spin-independent splittings in phys- 
ical units to be independent of the lattice spacing, indicating the absence of 
systematic errors from discretisation effects. Spin-dependent splittings are 
more sensitive to the lattice spacing and higher order corrections to the ac- 
tion; we discuss the size of these effects and what can be done to arrive at a 
physical result. 
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I. INTRODUCTION 



Accurate calculations of the hadron spectrum in lattice QCD require control of systematic 
errors. This has become a very important issue now that statistical errors have been reduced 
in recent years to the point where systematic errors can dominate the results reported. 

A major source of systematic error is that arising from the use of a space-time lattice 
with finite lattice spacing. All operators in the continuum Lagrangian must be replaced with 
discrete versions and discretisation errors consequently appear. This means that physical 
results (for example a mass in GeV) depend on the value of the lattice spacing. This is 
obviously wrong. Since the lattice is simply a regulator for the theory, physical results must 
not depend upon its value. One approach has been to extrapolate to zero lattice spacing 
for a 'continuum' result. This is difficult numerically, especially if the variation with lattice 
spacing is severe. However, recent progress in understanding discretisation errors and how 
to formulate an improved action has meant that we can obtain essentially continuum 
results at finite values of the lattice spacing. The lattice spacing dependence is reduced to 
such a low level that extrapolation is unnecessary. 

The spectrum of bottomonium bound states is one of the most accurate calculations 
that can be done on the lattice J3|. Since the b quarks are non-relativistic in these systems 
(y 2 /c 2 ps 0.1) a non-relativistic action can be used 0,0]. This allows a b quark propagator 
to be calculated on one sweep through the gluon field configuration with low computational 
cost. Multiple sources can be used on a single configuration because the bound states are 
much smaller than the volume of a typical lattice. Also sources for both ground and excited 
states can be used, allowing multi-exponential fits to hadron correlators and improving the 
confidence in the fitted masses. These techniques mean that very small statistical errors can 
be obtained and the improvement of systematic errors becomes a priority. In this paper we 
discuss the issue of discretisation errors for the bottomonium spectrum ||. 

The approach that we use, Non-relativistic QCD (NRQCD) is an effective field 

theory. Its Lagrangian is suited to a description of non-relativistic quarks since operators 
are classified according to the powers of v 2 /c 2 they contain, where v is the velocity of the 
heavy quark. The number of operators to be included can then be truncated at a fixed order 
in v 2 /c 2 and this is clearly a sensible thing to do if v 2 /c 2 <C 1. The renormalisability of QCD 
is lost in this process but physical results are still obtained by putting an explicit momentum 
cut-off into NRQCD. This cut-off should exclude relativistic momenta and thus be of the 
same order, or smaller, than the heavy quark mass. On the lattice this cut-off is provided 
by the lattice spacing, with Ma ~ 1. The excluded momenta cause renormalisation of the 
coefficients of the NRQCD operators when, say, lattice NRQCD is matched to full continuum 
QCD. The coefficients will be well-behaved and essentially cut-off independent provided that 
the cut-off is not too large. Any attempt to take the cut-off to infinity (lattice spacing to 
zero) will cause them to diverge as the non-renormalisability of the theory becomes apparent. 
Thus, no continuum extrapolation can be done for lattice NRQCD. However, as discussed 
above, a continuum extrapolation is not necessary for a suitably improved action. All that 
is necessary is to demonstrate lattice spacing independence of physical results. For the 
bottomonium spectrum from NRQCD this should be possible in a region of lattice spacing 
M^a ~ 1, and this is what we show in this paper. 

The size of discretisation errors will vary from one quantity to another. In general it is 
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to be expected that the coefficient of the dependence on a should represent some typical 
momentum scale appropriate to that quantity 0. For the light hadron spectrum this would 
then be a few hundred MeV. For heavy hadrons the scale of discretisation errors is likely to 
be larger. The scale is not set by the heavy quark mass since this is an irrelevant scale to 
the dynamics of the bound states. It is set rather by typical momenta exchanged inside the 
hadrons. For bottomonium these momenta are of order 1 GeV and so discretisation errors 
might be expected to present a problem on coarse lattices if the action is not improved. Here 
we report results with leading order (a 2 ) discretisation errors removed from spin-independent 
terms, but not from the spin-dependent terms (which are of lower order in the non-relativistic 
expansion) . 

Section 2 describes the lattice calculations and results at three different values of the 
lattice spacing. Section 3 discusses the scaling behaviour of spin-independent and spin- 
dependent splittings. Section 4 contains our conclusions. 



II. NRQCD CALCULATIONS AND RESULTS 

Quark propagators in lattice NRQCD are determined, in a single pass through the gauge- 
field configuration, from evolution equations that specify the propagator for t > in terms 
of its value at t = 0. We use here: 

aH \ n x / aH x " 



Hq is the kinetic energy operator, the lowest order (in v 2 /c 2 ) term in the Hamiltonian: 

A ( 2 ) 

Ho = - 2^. (2) 

The correction terms to the Hamiltonian that we include in SH are (9(w 4 /c 4 ). They comprise 
relativistic corrections to the spin-independent Hq as well as the first spin-dependent terms 
that give rise to spin-splittings in the spectrum. 



(A< 2 >) 
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The last two terms in SH come from finite lattice spacing corrections to the lattice laplacian 
and the lattice time derivative respectively ||. A is the symmetric lattice derivative and 
is a lattice version of the continuum operator J^Df. We used the standard traceless 
cloverleaf operators for the chromo-electric and magnetic fields, g>E and gB. The parameter n 
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is introduced to remove instabilities in the heavy quark propagator caused by the highest 
momentum modes of the theory. 

We tadpole-improve |J our lattice action by dividing all the gauge fields, U, that appear 
in E, B, and the covariant lattice derivatives fields by Uqp, the fourth root of the plaquette. 
This is most easily done as the U^s are read by the code that evolves propagators. Tadpole- 
improvement of the action allows us to work with tree-level values for the q's in 5H (i.e. 1) 
without, we believe, having to worry about large renormalizations [|TUJ. Hence our lattice 
action depends only on two parameters, the bare mass M 6 ° and the QCD coupling constant, g. 

Table | shows the parameters used in the calculations at 3 different values of (3. The 
configurations were all generated using the standard unimproved Wilson plaquette action 
and generously made available to us by the UKQCD collaboration JTlJ and by Kogut et al 



T2|| . The results described here at (3 = 6.0 agree with our previous results || but generally 
have higher precision, because of an increased number of sources on different time slices and 
the increased length of the lattice in the time direction. 

Once the quark propagators have been calculated it is straightforward to obtain anti- 
quark propagators and meson correlation functions. We used the standard interpolating 
operators described in || with source and sink 'smearing functions'. We worked in Coulomb 
gauge and took wavefunctions for smearing functions, either from a Richardson potential (/3 
= 6.0 and 6.2) or from a Coulomb potential (with modifications, (3 = 5.7). We took a ground 
state wavefunction and 2(1) radial excitations for S states at (3 — 6.0 and 6.2 ((3 = 5.7). 
For P states we used a ground state wavefunction and 1(0) radial excitations. We also used 
local sources and sinks which were delta functions for S states and combinations of delta 
functions for higher orbital excitations. In addition we looked at S-state mesons with small 
non-zero momenta. In the following discussion correlation functions at zero momentum will 
be denoted { 2S+1 Lj) a f 1 where a is the source smearing function and b the sink smearing 
function with I for a delta function (local operator), 1 for the ground state, 2 for the first 
excited state and so on. 

At j3 = 5.7 we summed over both initial quark spins. At (3 = 6.0 and 6.2 we saved CPU 
time by fixing the initial quark spin to +1/2, since the spin-flip operators in the Hamiltonian 
are suppressed with respect to H . We then used the strong correlations between different 
polarisations to obtain reduced errors on the spin splittings for P states [[|. 

We used 8 different spatial origins for our quark propagators at (2)(4)(1) different time 
slices at (3 = (5.7) (6.0) (6.2) to improve statistics. At (3 = 6.2 all 8 spatial origins were 
handled simultaneously by using Z(2) noise at each origin, 1 set per configuration [O . 

As described in || we used multi-exponential fits to the multiple correlation functions 
obtained by different combinations of source and sink. This allowed us to obtain ground 
state energies and one or two excited state energies. Two different types of fit were employed; 
the 'matrix' fit and the 'row' fit. The matrix fit used the matrix of correlators obtained with 
ground and excited state sources and sinks. The row fit used the row of correlators with 
ground and excited state sources and local operator sinks. We found the correlators with 
both local source and sink to be of very little use in fitting. We had a large number of 
measurements in every case and so did not run into problems with our covariance matrix, 
even for multi-exponential fits with several parameters. The two different fits gave consistent 
results within the errors that we quote. 

Final fitted values were chosen by monitoring the quality of the fits (Q) for given ranges 
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of fitting time, t m i n to t maX} as well as the stability of the fitted parameters. For a given 
type of fit, Q generally increases sharply with t min until it reaches a plateau. The first fit 
for which this happens is taken as the preferred value. In general the nth excited energy is 
taken as reliable from a fit to n + 1 exponentials. Table [II] shows the quality of our fitted 
results for a 2 exponential fit to 2 correlators for the 3 Si; ( 3 Si)ii and ( 3 Si)2z, at each (3 value. 
Notice how Q increases from small values of £ m j n as contamination from a third state dies 
away. Notice also how stable the fitted ground state energies are for a very large range of 
tmin values. 

It is interesting to study how the noise in the meson correlators changes with (3. We 
expect the ground state meson correlators ( 1 5'o) to have noise governed by the same mass 



as the signal since it is the lightest mass available [14]]. This means that the errors in 
an effective mass plot will not grow with lattice time. Figure 1 shows this clearly for the 
effective masses from the { 1 Sq)u correlator. The size of the errors at different values of (3 
reflects partly the different statistics available for the different sets of configurations (see 
Table [[]). If we multiply the errors at (3 of 5.7 by y/2 and at 6.0 by 2, for the different 
number of time origins (assuming these are independent), then the errors are in the ordering 
5.7 > 6.0 > 6.2. The error depends on the overlap between the squared correlation function 



and two 1 So particles |TJ]]. On coarser lattices, the local sink will provide a better overlap 
with two 1 Sq particles than on finer lattices and so we would expect the error to be larger. 
If instead we compare correlation functions in which the ground state smearing is applied at 
both source and sink, ( 1 5'o)ii, then the errors at all three values of j3 are very similar when 
adjustments for statistics are made as above. This reflects the fact that the noise should not 
change if the physical overlap with two 1 S , states doesn't change. 

For higher states than the ground state the noise grows exponentially with time according 
to the splitting between that state and the 1 5'o. Figure |2| shows this effect for the 1 P± 
correlator with ground state smearing at the source and a local sink. We expect the doubling 
time for the error to be In(2)/(1P — IS) = 1.6 GeV -1 , and this is roughly true at all three 
values of (3. Again the absolute size of the error at fixed physical time is very similar 
between all 3 (3 values when adjustments for the different statistics are made as above. 
However, on the coarser lattice many fewer lattice time points occur before the noise grows 
overwhelmingly large. 

Table [TTT| shows energies in dimensionless units obtained from our fits at each value of (3. 
3 exponential fits were used in general so the values for the 3S states should be used with 
some caution. 3 P is the spin average of the 3 Po,i,2 states defined by 

3 = 5M( 3 P 2 ) + 3M( 3 P) + M( 3 P ) 

9 ' 1 ' 

This is obtained by measuring spin splittings (see below) between the 3 P states and the 1 P\, 
and the energy of the l P\. 

Because the quark mass term is missing from our Hamiltonian the zero of energy becomes 



shifted so that the energies measured in the simulation and given in Table |TJ cannot be 
directly converted to hadron masses. Differences in energy can be converted directly to 
physical units using a value for the lattice spacing, but to obtain absolute masses we need to 
know the energy shift. It is sufficient to calculate an absolute mass for one meson only, and 
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the one for which the most accurate calculation can be done and compared to experiment 
is the l 3 ^x, the T. 

To calculate the absolute mass of the T we measure the dispersion relation from the en- 
ergy of meson correlation functions at small, non-zero momenta, and fit to a non-relativistic 
energy-momentum form: 

aE T (p) = aE NRiT + ^--C 1 ^^. (5) 



Enr is the energy at zero momentum normally measured (and given in Table |T|). C± is 
a constant to be obtained from the fit. aE{p) — aE NR is obtained accurately by a single 
exponential fit to the bootstrapped ratio of correlators at finite and zero momentum. We 
use the lowest 1 or 2 non-zero momenta in the fit. Table |T| shows the kinetic masses in 
lattice units, aMx, obtained at the 3 different values for f3 for the bare quark masses given 
in Table |. The value at (3 = 6.0 is taken from Ref. Jl and has not been recalculated on the 
Kogut et al configurations. 

P-wave spin splittings can also be obtained most accurately from ratio fits. Single expo- 
nential fits are performed to the ratio of appropriate polarisation components to maximise 
the correlations, as discussed above. Table [TTT] shows the results for different splittings in 



lattice units at the 3 different values of (3. Table [TV] gives a more detailed breakdown for 
different polarisation components of the IP fine structure at /3 — 6.0 and 6.2. For a given 
splitting there is no significant difference between different polarisation components so we 
average to get a final value and allow for variations in the error. We see no significant 
difference between 3 P2T states and 3 P2E states. This has been checked explicitly by taking 
ratios of those correlators. For the hyperfine splitting (T — rjb) we are able to extract both 
ground and excited splittings (for the first time in a lattice calculation) from simultaneous 
fits to the 3x3 matrix of correlators for the 1 Sq and 3 S\. These results are also given in 
Table fTJ. 

The wavefunction at the origin is calculated from the ratio of amplitudes of row and 
matrix fits as described in ref. 101. 



III. DISCUSSION 
A. Setting the scale 

One of the useful features of the spectrum of heavy quark bound states is that the 
splittings between radial and orbital excitations, spin-averaged, are to a good approximation 
independent of quark mass in the region between bottom and charm. Since not all the 
bottomonium fine structure has been seen experimentally, this statement relies to some 
extent on estimates of the spin splittings that have not been measured. However, since spin 
splittings are very small for bottomonium systems (~ 10% of spin-averaged splittings), we 
still expect little quark mass dependence for radial and orbital splittings in the region of 
M b when non-spin-averaged splittings are used. This allows us to set the scale from lattice 
calculations, independently of the requirement to tune the bare lattice quark mass to get 
the right kinetic mass for the T. 
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In Table [V| we show values for the lattice spacing, obtained by fixing various radial and 
orbital splittings to experiment (spin-averaging where possible), at the 3 different values of 
(3. These lattice spacing determinations are very accurate ones for setting the scale in the 



determination of a s |TJ]. Notice that different splittings at a given value of j3 give slightly 
different values for a~ l . This is a feature of the quenched approximation which we return to 
below. First we describe how the determination of a -1 is done and how errors are assigned. 

Since the NRQCD action of eq. ^ is corrected through 0(a 2 ) for discretisation errors we 
would also like to remove other 0(a 2 ) errors that come from using gluon fields generated 
with the simple plaquette action. Fortunately, these errors can be corrected for after the 
calculation. Perturbatively the correction appears as a shift to energies and is related to the 
wavefunction at the origin. It can then be written in terms of the hyperfine splitting for s 



states (for p states the shift is zero). We use [|TJ 

a AM g = — (a M b f a AM hyp , (6) 

with Mb set to 5 GeV. The resulting shift to the splittings of Table |V| is given in column 4. 
This shift is added to the splittings in lattice units before they are divided into the physical 
splitting to obtain a -1 . For the shift for 2S states we use the ratio of 2S to IS hyperfines 



determined at (3 — 6.0 and given in Table [II. The statistical error in the splitting is then 



inflated by aAM g /2 before calculating the statistical error in a~ l given in column 5. 

This determination of a -1 also has systematic errors, as in all lattice determinations. 
We attempt here to quantify the errors relevant to our calculation. There are two sources, 
physical and unphysical. The errors from higher order relativistic corrections which have 
been ignored are physical and will give the same percentage error at all values of the lattice 
spacing. The unphysical systematic errors come from higher order discretisation corrections 
that have not been included. These are all much larger at the coarsest lattice spacing than 
elsewhere. 

As a guide to estimating the size of these errors we have estimated, using a potential 
model [16 1, the size of the shifts in energy caused by the relativistic and discretisation 



corrections of 5H that we have included and compared those to the results of NRQCD 
calculations 0. The potential model estimates for the sum of relativistic and discretisation 
corrections in 5H at j3 — 6.0 give a resulting shift to the bb IP — IS splitting of —10 MeV, 
made up of —10 MeV from the relativistic corrections (10 MeV from the p 4 term and —20 
MeV from the D-E term) and cancelling contributions each of around 10 MeV in magnitude 
from the two discretisation corrections ||16|| . The 25 — 15 splitting has smaller shifts because 
the expectation values on which the shifts depend are more similar for the 2S and 15* than 
the IP and 15*. The expected result from adding all the terms is still —10 MeV. The NRQCD 
results show a IP — IS splitting that is 20(30) MeV larger without SH than with, and a 
2S — IS splitting that is 5(25) MeV larger. This is in good agreement with expectations, 
albeit with statistical errors that are too large to show a clear effect. However, if the shifts 
had been much larger than the estimates, they would have been visible above the noise. 

From this result we can extrapolate to the size of relativistic and discretisation corrections 
that we have not included. The terms that we can estimate most readily are those involving 
powers of quark momenta since these are easy to relate to lower order terms. The terms 
involving chromo-electric and magnetic fields and those terms with a structure that appears 
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for the first time at higher order are much harder and we have not estimated these. There 
is no reason, however, to suppose from our study above (which compares the p A and Darwin 
terms), that these terms should be any larger than the ones we can easily estimate. 

Higher order relativistic corrections would appear as v 6 terms in 6H, i.e. 0(v 4 ) relative to 
the leading terms. The percentage error we expect is then naively 10% 2 = 1%. On the other 
hand, the estimates above using potential models |16| , of the v 4 spin- independent relativistic 
corrections that we have included show them each to be less than half of the 10% (= 50 
MeV) naively expected. In addition we actually need the difference between the corrections 
for, say IP and IS, to get the shift in the splitting. This indicates that higher order 
corrections could be smaller than 1% too. Another type of similar higher order correction 
is that from radiative corrections to the q coefficients beyond tadpole-improvement. These 
should appear at the level of a s at an ultra-violet scale times the v 4 relativistic corrections, 
giving 0.5% at (3 — 5.7 and less at higher (3 values. To encompass both these higher order 
physical corrections a 1% error is given as the second error in column 5 of Table [V]. 

The NRQCD action, eq. ^|, includes the leading a and a 2 corrections which appear in 
power counting form as p 2 a 2 v 2 and Kav 2 relative to the leading order v 2 term (Hq). p and 
K are a typical momentum and kinetic energy associated with the bound state. Estimates 
of these terms using potential models |TB| yield shifts of 40 MeV at (3 = 5.7, 10 MeV at (3 = 
6.0 and 5 MeV at (3 = 6.2, in the IP — IS splitting. Higher order discretisation errors not 
included could be radiative corrections to those included (beyond the tadpole-improvement 
of these terms which has been done) and we can estimate these as a s (7c/a) [|TIJ times the 
leading errors. This gives 7 MeV at j3 = 5.7, 1.5 MeV at (3 = 6.0 and 0.7 MeV at /3 = 6.2. 
Higher order terms in a (such as p 4 a 4 v 2 terms) would give a percentage effect roughly the 
square of the leading order terms, i.e 4 MeV at (3 = 5.7, 0.2 MeV at (3 = 6.0 and essentially 
at f3 = 6.2. We should also consider the first discretisation corrections to the first relativistic 
corrections i.e. terms of order p 2 a 2 v A . These are very similar looking terms to the p 4 a v 2 
terms and so we can use this to estimate their size. Including powers of Ma and numerical 
factors we get a similar size correction at /3 — 5.7. At higher values of (3 these corrections 
are actually more important than those which are higher order in a but not suppressed by 
powers of v 2 . However, all discretisation corrections become smaller at higher (3 so they are 
still negligible; 0.5 MeV at (3 = 6.0 and at (3 — 6.2. The third error given in column 5 of 
Table |V] is then conservatively estimated by the sum of the three errors given above from 
discretisation errors in the NRQCD action. It is interesting to note that the discretisation 
corrections are smaller than the anticipated higher order relativistic corrections except at (3 
= 5.7. Note also that the statistical errors are generally larger than the systematic - to see 
any effect from including higher order terms we would have to improve our statistical error 
significantly. Higher order discretisation corrections from the gluon action are anticipated 
to be negligible given the size of the 0(a 2 ) correction in Table [V|. 

The a -1 determination at (3 = 6.0 agrees with our previous determination B and at (3 



6.2 agrees with previous UKQCD results [17 
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B. Determining the quark mass 



As discussed above, experimental spin-averaged radial and orbital splittings are very- 
insensitive to the value of the quark mass. However, in the quenched approximation there 
is some mass dependence for these splittings on the lattice |1§ , and this is increased if spin- 



averaging is not done. It is therefore true that a tuned bare quark mass is necessary to get 
the right radial and orbital splittings. The spin splittings are much more sensitive to the 
quark mass (roughly as its inverse), and it is essential to tune the quark mass to get these 
correct. The way in which we tune the quark mass is to adjust it until the kinetic mass of 
the T agrees with experiment. Table |VT| shows the kinetic mass of the T (given in lattice 
units in Table |T|) in GeV, for each of the 3 values of (3, using a~ l from Table |V|. There is 



reasonable agreement with the experimental result 9.46 GeV in each case, provided that the 
value for a~ l is taken from the T — T splitting. The systematic errors in the determination 
of aMki n are at the 1% level from the same sources as systematic errors in the determination 
of a~ l . This is smaller in every case than the statistical error, dominated by the statistical 
uncertainty in a -1 . 

The difference between the kinetic mass of a meson and its energy at zero momentum is 
calculable in perturbation theory ||19|| . The formula which relates E^r and the kinetic mass 
M is : 

M = 2(Z m M° - E ) + E NR (7) 

where Z m is the mass renormalisation and E Q the energy shift. Table |VII| gives values for Z m 
and E appropriate to the different bare masses used at each value of (3. Tadpole-improved 
lattice perturbation theory has been used for these parameters and the scale of a s set using 
the BLM scheme ||19| . The values obtained for My from the perturbative expression, (0), 



are given in lattice units in the sixth column and should be compared with the results for 



aMkin in Table [VI]. There should be agreement at all values of (3 independent of whether 
the quark mass is well tuned to that appropriate to the b or not, and we see that there is. 
The perturbative error is at 0(a 2 s ) and is taken here as the square of the 0(a s ) term. Note 
that the relationship in eq. [7] is well defined perturbatively. Non-perturbative values for the 
shift between Enr and M can also be measured on the lattice for, say, T and used, divided 



by 2, for B physics |2C 



We can use the information above to determine the b quark mass [21] in two independent 



ways and ask whether the renormalised b quark mass we obtain scales from one value of (3 



to the next. These results will be presented elsewhere [22| 



C. Radial and Orbital splittings 



From the results of Table [TTJ we can investigate the scaling of spin-averaged splittings 
across the range of (3 values we have used. Our analysis of systematic errors in section 3.1 
already implied that we do not expect to see violations of scaling in these quantities. 

Figure |] shows with circles the ratio of our xE~ T splitting to the p mass, obtained by the 
UKQCD collaboration at the same three values of f3 []TT| , p3| . The p mass was calculated using 
an action with 0(a) errors removed using a clover term with tadpole- improved coefficient; i.e. 
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a similar philosophy to that used here to remove discretisation errors from the T splittings. 
The remaining errors in the p determination may then be of 0(a 2 ). For the IP — IS splitting 
we have removed the 0(a 2 ) gluonic errors so, as discussed earlier, the remaining errors are 
higher order in a and/or suppressed by powers of v 2 . The plot is very flat showing that good 
scaling is obtained. We can rule out scaling violations in this ratio with a scale larger than 
a few hundred MeV, and it would be very unlikely for the p mass and the T system to have 
identical and cancelling large scaling violations. The plot also shows in the comparison to 
the GF11 results, that if an unimproved calculation is done of the p mass [24], an absence 
of scaling is quite evident (in this case linear in a) . 

The ratio A(xE — Y)/m p shows a big discrepancy with experiment. This we believe is an 
error from the quenched approximation. The scales intrinsic to a light hadron system and 
the T are quite different and the coupling constant does not run correctly between these 
scales in the quenched approximation. Again, these errors are masked by discretisation 
errors for an unimproved light hadron spectrum. 

Figure § also shows the ratio of the \b — T splitting to the A parameter of QCD in the 
V scheme ||. Again good scaling is seen over this limited range of lattice spacings. 

Figure |4] shows the dimensionless ratios of various splittings within the T system (spin- 
averaged where possible) as a function of lattice spacing. The ratios shown are for the 
T — T, T' — T and the h b — T (experiment for the latter case using \ b for h b ). The 
ratios are constant as a function of lattice spacing, although the statistical errors are rather 
large for the 3S and 2P cases. This means that we can interpret our results as continuum 
results. Note that our current value for the energy of the 2 3 S\ at (3 — 6.0 has dropped by 
la from our previous calculation || with slightly smaller statistical errors. Our 1 1 P 1 energy 
has not changed and the error has fallen by a factor of 2. Our best value for the ratio 
(2S — IS) /(IP — IS) is now 1.35(6) and the discrepancy with experiment remains (1 — 2)a, 
controlled by the uncertainty in the 2 3 S' 1 energy. The discrepancy with experiment for the 
ratios involving the 3S and 2P could be affected by finite volume errors as well as quenching 
errors. 



D. Spin splittings 



Spin splittings in lattice units are given in Table M. This includes, for the first time, 
the radially excited hyperfine splittings, T — T] b and T — r\ b at (3 = 6.0 where our results 
are most precise. No hyperfine splitting has yet been seen experimentally for the T system, 
and a radially excited hyperfine splitting is still awaited for charmonioum. We find the 2S 
hyperfine splitting to be about half that of the ground state here. A proper analysis however, 
has to study the lattice spacing dependence and the effect of unquenching on this ratio. 

In this section we investigate the dependence in physical units of the ground state spin 



splittings on the lattice spacing. We also compare to recent UKQCD results |25j in which 
a systematic study at next order in relativistic and discretisation corrections has been done 
for spin-dependent terms, i.e. those which affect the spin splittings, only. By comparing the 
scaling of our results with their results we can untangle to some extent the difference between 
the discretisation errors at fixed a, which are unphysical, and the relativistic corrections, 
which are physical. Whether or not the spin-dependent relativistic corrections (i.e. terms of 
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OiMv 6 )) are sizeable or not is important since these are terms outside the scope of potential 
models. 

To investigate scaling we show in Figures |5] and ^ various spin splittings in MeV setting 
the scale from the T — T splitting as a function of a 2 . The T — T splitting is used since this 
is the one for which our quark masses are closest to being tuned and this is important for spin 
splittings [ 26|j . The results from Table |T| are given as squares. We include in the Figures 



the direct uncertainty in the spin splittings from statistical errors in the lattice spacing. 
There is an additional error, which we have not included, from the uncertainty in the kinetic 
mass from the error in the lattice spacing. This uncertainty leads to an uncertainty in the 
tuned quark mass which affects the spin splittings. If the spin splittings vary as the inverse 
quark mass (roughly true at least for the hyperfme splitting), then this second uncertainty 
is actually equal in size, correlated with and in the same direction as the error that we have 
included. This would lead to an approximate doubling of the error bars in the Figures. 

Clear dependence on the lattice spacing is visible for the hyperfme splitting, T — 
For the IP fine structure it is less clear because statistical errors, particularly at (3 = 6.2, 
are so much higher. Lattice spacing dependence in our results is not surprising since the 
Hamiltonian we have used has only leading order spin-dependent terms with no corrections 
for discretisation errors in these terms. These errors are expected to be 0(a 2 ) relative 
to the leading term and therefore in relative terms significantly larger than for our spin- 
independent splittings. In addition, in a potential model picture the fine structure is provided 
by potentials which are generally of much shorter range than the central potential so we 
would expect the fine structure to have a harder scale for discretisation errors than the 
radial and orbital splittings. If we parameterize the scaling violations in the hyperfme 
splitting by 0: 



splitting I splitting\ 
2S-1S = \2S-1S j. 



1-//V), (8) 



we find the parameter p, to take a value around 0.8 GeV for the squares. 

For the IP fine structure we have the advantage that we can compare to experimental 
results as well as comparing results at different values of the lattice spacing. In Figure |6| we 
notice a clear disagreement with experiment for the squares in that both the overall scale of 
splittings is too small and the ratio 



is too large. Similar results were noticed in our charmonium spectrum calculation p7| . The 
disagreement with experiment over p was suggested there as a discretisation effect and results 
here tend to confirm that discretisation errors do have an effect in this ratio. Our value for 
p at beta=5.7 is 1.4(3) and at 6.0, 1.1(3). The experimental value is 0.66. The comparison 
between the squares at (3 = 5.7 and 6.0 shows that scaling violations are just about visible for 
the 3 Po and 3 P2 above the statistical errors. The slope of the scaling violations is similar to 
that for the hyperfme splitting (but with large errors). Also in agreement with the hyperfme 
splitting we can see that adding discretisation corrections would increase the splittings and 
this would make the IP fine structure closer to experiment. 
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On Figures [5| and ^| we also show results from Manke et al [25 1 at /3 = 6.0 and 5.7 with 
diamonds. They have used the same bare quark masses as us at both values of (3 but included 
higher order spin- dependent relativistic (0(v 6 )) and discretisation {0(a 2 v )) corrections to 
the action ||, all tadpole-improved with ti p. We take their spin splittings in lattice units 
and convert to physical units using our T — T splitting since, as explained earlier, this is 
the a -1 for which the quark masses are best tuned. Our result for this splitting should be 
the same as theirs since they have added only spin-dependent terms |2J| and ours is more 
precise. The lattice spacing dependence of their results should be much reduced over ours, 
and we see that it is. By comparing their results and ours, we can see that the relativistic 
corrections act in the opposite direction to the discretisation corrections, in agreement with 
the findings of |29[ . 

By comparing the results of [25J with ours we can attempt to untangle the size of v 6 
relativistic corrections to the fine structure. First, however, we should also include radiative 
corrections to the leading order terms, beyond tadpole improvement. These are of 0(a s v 4 ) 
and therefore of a similar size to the relativistic v 6 corrections. The Cj in equation 3 take 
the form 1 + c\ a s (q*) + . . . with of* of 0(1). Preliminary perturbative results are available 
p0| for c 4 , for the a ■ B term in the action, and confirm this expectation. In perturbation 
theory the hyperfine splitting then picks up a factor of the square of this correction, so it 
is quite sensitive to c 4 . Figure [5] shows with fancy squares our results for the hyperfine 
splitting rescaled by (1 + a;p(3.4/a)c 4 ) 2 . Values for ap were taken from [|l5j. C4 increases 



as Ma increases, which gives a larger correction at low (3. Since this is also the direction in 
which ctp(3A/a) increases the correction at (3 = 5.7 is large, 50%. At (3 — 6.2 the correction 
from c 4 is 20% [3lf[. The fancy diamonds show the rescaled results from |25| in which we 



have simply rescaled the leading order piece i.e. we have applied the same shift in MeV to 
these results as to ours. It is clear from this Figure that the radiative corrections act in 
the same direction as discretisation corrections discussed earlier and absorb some of them, 
softening the value of \x in equation [||. 

For the hyperfine splitting, we will assume that the discretisation corrections at (3 = 6.0 
in the (v 6 , a 2 v 4 ) action J25J are of 0(a 4 ) and therefore negligible already at (3 = 6.0. For 
the (f 4 ) action (this work) we will take the result at (3 — 6.2 as being the least affected 
by discretisation errors (we don't expect more than a 5% effect here). The difference then 
between the result in physical units between (3 = 6.2 (v A ) and f3 = 6.0 (v 6 , a 2 v 4 ) represents 
the effect of relativistic corrections. From Figure |5j with fancy squares and diamonds this 
gives a 15% decrease in the hyperfine splitting and from the plain squares and diamonds a 
20% decrease, consistent with the expectation of around 10% from power counting in v 2 /c 2 . 

For the p fine structure, the radiative corrections are not available and they would in 
any case be much harder to apply without doing the full lattice calculation including them. 
We therefore just compare the results with (diamonds) and without (squares) the v 6 and 
a 2 v 4 corrections in Figure ||. The lattice spacing dependence is apparently reduced by the 
a 2 v A corrections, although the statistical errors make it hard to quantify this result. The 
relativistic corrections again act to reduce the splittings, now to something like 50% of 
experiment. The size of the relativistic corrections is significant, possibly several times the 
naive estimate of 10%. They could be offset by radiative corrections, however, in a fully 
consistent calculation that includes the next order in v 2 , a 2 and a over our calculation. 
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Radiative corrections which increase the spin splittings, as for the hyperfine case, would 
improve the agreement with experiment, but unless these are surprisingly large, the splittings 
will remain too small in the quenched approximation. Initial unquenched results for the v 4 
action and the v 6 , a 2 v 4 action |32| have indicated an increase in the p spin splittings, but 
currently statistical errors are too large to clearly determine the size of unquenching effects. 

Ref. |25j points out that the (v 6 , a 2 v 4 ) action has a much improved value of p (eq. H) 
over the (t> 4 ) action. They quote 0.56(19) at /3 — 6.0 in the quenched approximation, now 
in agreement with experiment, confirming the discretisation effects for this ratio discussed 
above. In fact, we do not expect this ratio to agree exactly with experiment in the quenched 
approximation. An analysis of p for the simple Cornell potential |33] shows that for a 
quenched a s which is too small, the value of p is reduced. With better statistics it may 
become clear that the quenched value for p is actually smaller than the experimental result. 

Another way to study the sensitivity of the fine structure to higher order terms in the 
action is to change the value of the tadpole coefficient, u . Several calculations ]T| have 
indicated that the average link in Landau gauge, m l is to be preferred over the fourth root 
of the plaquette, wop, that we have used here. Changing uq is equivalent to resumming 
various terms to all orders in the Cj. It is again easy to rescale the hyperfine splitting after 
the simulation to account for a different value of u , by counting the powers of u that appear 
in front of the a ■ B term. The hyperfine splitting is very sensitive to uq, scaling as 1/u®, 
when allowing for uq factors in H which change the renormalised quark mass. The expected 
scaling is confirmed by calculations at different values of uq || (although here a fixed bare 
quark mass was used). The hyperfine splitting increases most on using uql at (3 = 5.7. There 
the ratio of uqp/uql is 1.045; at (3 — 6.0 it is 1.020 and at (3 — 6.2, 1.013 f3"4]| . The hyperfine 



splitting with m l then has somewhat reduced lattice spacing dependence. Note that c 4 will 
be different for the action with m l, since a different amount of the radiative correction has 
been absorbed by the tadpole factors. C4 in fact is larger, by 0.24, being the cube of the 
ratio of the uq factors to 0(a). The size of relativistic and discretisation corrections will 
change by differing amounts on changing uq because they have different effective powers of 
Uq in their coefficients. 

From the above we can conclude that systematic errors in the fine structure still exceed 
10%, even in the quenched approximation. To improve the accuracy the v 6 and a 2 v 4 action 
must be used, with radiative corrections included for the leading spin- dependent terms. For 
the hyperfine splitting the sensitivity to the q is such that the O (a 2 ) term also needs to be 
known (or at least bounded). The determination of the q is in progress []30]]. More accurate 
meson kinetic masses and lattice spacing determinations are also required so that the quark 
mass can be tuned more accurately, since the spin splittings are sensitive to this. 

The wavefunction at the origin is another short distance quantity but one which does 
not vanish as the spin-dependent terms are switched off. In Figure |7| we show results for 
ip(0) for the T and T as a function of a 2 , again setting the scale from the T — T splitting 
since the wave function is sensitive to quark mass. The lattice spacing dependence does 
not seem severe but the statistical errors are large. They are dominated by those from a -1 
which here is raised to the 3/2 power. The bursts show experimental values based on the 
leptonic widths and the naive Van Royen-Weisskopf formula, 
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v,- = W4%^, (io) 



where M v is the vector meson mass. There will be radiative corrections to this formula of 
0(a s ). For potential model wavefunctions these radiative corrections have coefficient 16/37T 
but this does not have to be the same for NRQCD. Calculations of the renormalisation 



factors have only been done for lower order actions than the one used here |35[]. There will 
also be terms of 0(1/Mb,a s /Mb) from mixing of the vector current with higher dimension 
operators (in a similar way to the mixing for J'b which has been worked out in some detail 

0)- 



A lot of the radiative corrections cancel in the ratio of leptonic widths for the IS* and 2S 
states. Experimentally this ratio is 2.5(2). Our result at (3 = 6.0 simply from the ratio of 
^(O)! 2 is 1.7(4). We expect errors in the quenched approximation which would suppress the 
IS wavefunction at the origin relative to the 2S, so that it is not surprising that our result 
is somewhat smaller than experiment. We cannot yet say whether the reduction observed is 
reasonable. 

To obtain results for the real world from simulations in the quenched approximation 
and at 2 flavours of dynamical quarks, extrapolations in Nf (= the number of flavours of 
dynamical quarks) will be required. For this it is important to have scaling quantities, as 
we have discussed here. It is also necessary to use quantities which are sensitive to the 
same number of dynamical fermions. This means that for a given splitting we should set 
the scale using a quantity which has similar internal momentum scales. For example, for 
the fine structure it may not be optimal to use the softer radial and orbital splittings to set 
the scale. Typical momentum scales for radial and orbital splittings are around 1 GeV so 
they can be expected to 'see' 3 flavours of dynamical fermions in the real world ||15|| . The 
fine structure may instead 'see' 4 flavors and this would imply that an extrapolation in Nf 
can only be done in terms of a quantity which also sees 4 flavors. 

In Figure [5] we show a plot of the hyperfine splitting in MeV against the square of the 
lattice spacing using the 1 3 P 2 — 1 3 P splitting to set the scale. Statistical errors are now 
very large, despite using the most accurate of the P fine structure splittings. Much higher 
accuracy will be needed to use this ratio for Nf extrapolations. 

IV. CONCLUSIONS 

We have calculated the spectrum of bottomonium bound states at 3 different values of the 
lattice spacing using NRQCD in the quenched approximation. Ratios of spin-independent 
splittings are constant in this region, independent of the lattice spacing. This is a necessary 
requirement for the results to make physical sense. There is no necessity to extrapolate 
to vanishing lattice spacing (impossible anyway for NRQCD). The constancy of the results 
implies that they can now be used in extrapolations to real-world Nf values when combined 
with results on dynamical configurations |15| , |32| . 

For spin-dependent splittings there is some visible (unphysical) lattice spacing depen- 
dence. This is not surprising since our spin splittings are only accurate to leading order, 
given the terms in the Hamiltonian that we have used. A comparison to a calculation 
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which does include higher order relativistic and discretisation corrections to spin-dependent 
terms shows that the relativistic spin-dependent corrections could be as large as 20% for 
the hyperfine splitting where they can be quantified. They act to reduce spin-dependent 
splittings in all cases, and this increases the disagreement with experiment for the P fine 
structure in the quenched approximation. 

For extrapolations to real world values of Nf we may need to employ different techniques 
for fine structure to that for radial and orbital splittings since different momentum scales are 
involved and this may lead to a different number of dynamical flavors being 'seen'. Statistical 
errors are currently too high to use a P spin splitting to set the scale in such a programme. 
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FIGURES 



FIG. 1. Effective mass plots for the ( 1 Sq)u correlator at all three values of /3, in order with (3 
= 5.7 at the top. The time axis has been converted to physical units of GeV -1 using the xE — T 
splitting to set the scale (Table |V|). 

FIG. 2. Effective mass plots for the ( 1 -Pi)i; correlator at all three values of (3, in order with (3 
= 5.7 at the top. The time axis has been converted to physical units of GeV" 1 using the xE — T 
splitting to set the scale (Table |v|). 

FIG. 3. Dimensionless ratios of the xt — T splitting to the parameter Ay (diamond) and to 
the UKQCD p mass [|ll]] (circle) against the lattice spacing in fm (set by the x& — T splitting). 
Results using the GF11 p mass |24| are given by squares. The burst represents the experimental 
value for A(xE — T)/m p . 

FIG. 4. Dimensionless ratios of various splittings to the Xb — T splitting against the lattice 
spacing in fm (set by the Xfe — T splitting). Circles represent the ratio for the T" — T splitting 
(experiment short dashes) and crosses for the h b — T (experiment using x'b for h' b dot-dash). The 
diamonds show the T — T ratio with a 2 gluonic corrections (as described in the text) and the 
squares uncorrected results (experiment dashed line). The squares and crosses have been offset 
slightly in the horizontal direction for clarity. 



FIG. 5. The hyperfine splitting in MeV using the T — T splitting to set the scale, vs a 2 in fm 2 . 
Plain squares indicate our results from Table III. The diamonds indicate the results from ref. [25] 
using a higher order action. They are at matching values of (3 but offset slightly for clarity. Fancy 
squares indicate our results, rescaled by the square of C4 calculated to 0{a) in |3(|. The fancy 
diamonds indicate the result from ref. |2f| shifted by the same amount as our results to account 
for radiative corrections to the a ■ B term. The error bars shown are statistical only and include 
some of the error from the uncertainty in the scale (see text). The x axis errors from uncertainty 
in the scale are shown only for the squares for clarity. 



FIG. 6. The splitting in MeV between various 1 3 P states and the 1 3 P using the T - T 
splitting to set the scale, vs a 2 in fm 2 . Squares indicate our results, diamonds the results of [25] 
at matching values of (3 but offset slightly for clarity. The bursts indicate the experimental values. 
The error bars shown are statistical only and include some of the error from the uncertainty in the 
scale (see text). The x axis errors are not shown for clarity - they are of the same size as those in 
Figure || 
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FIG. 7. The wavefunction at the origin in GeV 3 / 2 vs a 2 in fm 2 setting the scale from the 
T' — T splitting. Points for the l 3 S*i are shown as squares and 2 3 5i as diamonds. The errors 
shown are statistical only and the errors on the x axis are shown only for the squares. The bursts 
show experimental points derived from the leptonic widths of the T and T' using the naive Van 
Royen-Weisskopf formula. 

FIG. 8. The hyperfme splitting in MeV using the 1 3 P2 — l 3 -fb splitting to set the scale, vs a 2 
in fm 2 . Squares show our results, diamonds the results of |2E|. 



19 



TABLES 



TABLE I. The parameters used in calculations at 3 different values of the QCD coupling, 
(3 = 6/g 2 . 









n 


Uop 


V 


# confs. 


# sources 


collab. 


(3 = 


5.7 


3.15 


1 


0.861 


12 3 x 24 


200 


8x2 


UKQCD 


(3 = 


6.0 


1.71 


2 


0.878 


16 3 x 32 


149 


8x4 


Kogut et al 


13 = 


6.2 


1.22 


3 


0.885 


24 3 x 48 


216 


8 (Z(2)) 


UKQCD 



TABLE II. Fitted results for a two exponential fit to two 3 Si correlators, ( 3 <Si)i/ and ( 3 <Si)2Z- 
Fitted energies in lattice units are given with errors as well as the Q value for the fit at all three 
values of (3. 



(3 




5.7 






6.0 






6.2 




rnin 


E x a 


E 2 a 


Q 


Eia 


E 2 a 


Q 


Eia 


E 2 a 


Q 


2 














0.3128(3) 


0.517(3) 


0.00 


3 


0.5186(6) 


0.888(5) 


0.01 








0.3130(3) 


0.508(3) 


0.00 


4 


0.5186(6) 


0.901(7) 


0.37 








0.3131(3) 


0.504(4) 


0.00 


5 


0.5186(6) 


0.91(1) 


0.64 


0.4540(2) 


0.717(2) 


0.00 


0.3131(3) 


0.499(4) 


0.02 


6 


0.5188(6) 


0.93(1) 


0.73 


0.4539(2) 


0.710(3) 


0.14 


0.3132(3) 


0.499(5) 


0.04 


7 


0.5188(6) 


0.95(3) 


0.70 


0.4539(2) 


0.708(3) 


0.14 


0.3133(3) 


0.493(6) 


0.10 


8 


0.5187(6) 


0.93(4) 


0.65 


0.4539(2) 


0.705(4) 


0.13 


0.3133(3) 


0.491(7) 


0.08 


9 


0.5186(6) 


0.93(6) 


0.70 


0.4539(3) 


0.697(6) 


0.15 


0.3132(3) 


0.488(8) 


0.07 


10 








0.4539(3) 


0.697(6) 


0.14 


0.3131(3) 


0.49(1) 


0.07 


11 








0.4539(3) 


0.690(7) 


0.14 


0.3130(3) 


0.49(1) 


0.06 


12 








0.4538(3) 


0.68(1) 


0.13 


0.3130(4) 


0.48(1) 


0.05 


13 








0.4537(3) 


0.69(1) 


0.30 


0.3131(4) 


0.49(2) 


0.05 


14 








0.4537(3) 


0.70(2) 


0.24 


0.3133(4) 


0.50(2) 


0.04 
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TABLE III. Fit results for dimensionless bb energies and splittings, aE and aSE for the 
quenched approximation, at three different values of (3. Below these are given the kinetic mass (see 
text) and the wavefunction at the origin, in lattice units. 



3 


5.7 


6.0 


6.2 


Energies, aE : 








l 1 ^ (m) 

^ U \'IO/ 

l 3 5i (T) 


0.5029(5) 
0.5186(6) 


0.4415(3) 
0.4537(5) 


0.3028(2) 
0.3132(3) 


2% (%) 


0.92(3) 


0.678(8) 


0.478(6) 


2 3 5i (T'l 
3 3 5i (T") 


0.94(4) 


0.686(8) 
0.83(3) 


0.488(8) 
0.65(4) 


l 1 ^! (h h ) 
2 1 Pi (h' h ) 


0.843(6) 


0.627(3) 
0.823(14) 


0.438(5) 
0.60(7) 


1 3 P (w) 


0.845(6) 


0.628(3) 


0.440(5) 


Splittings, a&E : 








2 3 5"i - 2 1 S'o 
3 3 Si - 


0.01575(8) 


0.01237(14) 

0.006(1) 

0.005(3) 


0.01038(14) 


13 p _ l3 p 


020(2) 


0147(17) 


021(7) 


l 3 Po - l 3 Pi 


0.011(2) 


0.0078(15) 


0.010(6) 


1 3 P1 - 1 3 P 


0.0079(5) 


0.0069(12) 


0.010(7) 


l 1 ^! - 1 3 P;l 


0.003(2) 


0.0028(6) 


0.003(4) 


1 3 P 2 - 1 3 ^P 
1 3 P - l 3 Pi 


0.0059(7) 
0.0052(12) 


0.0042(5) 
0.0036(8) 


0.0058(24) 
0.0046(37) 


1 3 P - 1 3 P 


0.0137(11) 


0.0105(10) 


0.0148(52) 


1 3 P- l X Pi 


0.002(2) 


0.0008(8) 


0.0018(28) 


Kinetic mass, aM^in '■ 








M(l 3 5i) 


7.06(7) 


3.94(3) 


2.89(3) 


Wavefunction at the origin : 



|^(0) 
|^(0) 



| a 3 / 2 for l 3 Si 

|a 3 / 2 for 2 3 5i 

|a 3 / 2 for 3 3 5i 

|a 3 / 2 for 1% 

| a 3 / 2 for 2 1 S* 

| a 3 / 2 for 2 1 6* 



0.385(5) 
0.30(2) 



0.1525(14) 

0.118(14) 

0.19(3) 

0.1621(13) 

0.115(12) 

0.20(3) 



0.1116(12) 



0.1314(12) 
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TABLE IV. Individual IP spin splittings in lattice units with quantum numbers, polarizations 
and smearing combinations specified, along with the final value used in Table III. E and T stand 
for the different lattice representations of the continuum spin 2 operator. 



splitting 


(5 = 6.0 


result 


= 6.2 


result 


:i P 2 T yz {l,loc) 


- l P iy (l,loc) 


0.0060(7) 


0.0050(14) 


0.0072(47) 


0.0076(42) 


^2^(1,1) 


- x Piy(l,l) 


0.0048(7) 








3 P 2 T zx {l,loc) 


- 1 P 1 x(l,loc) 


0.0065(8) 




0.0060(52) 




3 P 2 T zx (l,l) 


- X PXX(1,1) 


0.0045(7) 








3 P 2 E zx (l,loc) 


- 1 P 1 z(l,loc) 


0.0050(8) 




0.0083(41) 




3 P 2 E ZX (1,1) 


- 1 P 1 z(l,l) 


0.0046(7) 








3 P 2 E yz (l,loc) 


- 1 P 1 z(l,loc) 


0.0048(8) 




0.0072(40) 




3 P 2 E yz (l,l) 


- X PIZ{1,1) 


0.0039(5) 








1 P 1 x(l,loc) - 


- 3 P iy (l,loc) 


0.0032(4) 


0.0028(6) 


0.0027(38) 


0.0028(43) 


l Pix{l,l) - 


- 3 ^iy(i,i) 


0.0026(5) 








l P iy (l,loc)- 


- 3 P 1 x(l,loc) 


0.0031(3) 




0.0039(42) 




l Piy{l,l)- 


- 3 Pix{hl) 


0.0024(4) 




0.0019(63) 




1 P 1 z(l,loc) - 


- 3 Po(l,loc) 


0.0097(10) 


0.0097(10) 


0.0139(49) 


0.013(5) 


^(1,1)- 


- 3 Po(l,l) 






0.0120(64) 





TABLE V. a -1 values at the three different values of (5 - the first error given is statistical, the 
second, systematic from higher order relativistic corrections and the third, systematic from higher 
order discretisation corrections. The experimental values for the splittings are 440 MeV (xE ~ T) 
and 563 MeV (T' - T). 

(3 splitting aAM aAM g cT 1 (GeV) 

5?7 Xb- T 0.326(6) -0.015 1.41(4)(2)(5) 

T'-T 0.42(4) -0.007 1.36(13) (2) (4) 

6.0 X6-T 0.174(3) -0.004 2.59(5)(3)(1) 

T'-T 0.232(8) -0.002 2.45(8)(3)(1) 

6.2 Xb-T 0.127(5) -0.002 3.52(14) (4) (0) 

T'-T 0.175(8) -0.001 3.24(15)(4)(0) 



TABLE VI. Values for the T mass at the three different values of j3, using different prescriptions 
for a -1 . The first error comes from the statistical error in aM^in, the second, the statistical error 
in a -1 from Table 0. The experimental value for the T mass is 9.46 GeV. 



p 


aM 6 ° 


aMkin 


M T (GeV) 

a (Xi,-T) 


M T (GeV) 


5.7 


3.15 


7.06(7) 


9.95(10)(30) 


9.60(10)(90) 


6.0 


1.71 


3.94(3) 


10.20(8)(20) 


9.65(7)(30) 


6.2 


1.22 


2.89(3) 


10.17(10)(40) 


9.30(10) (40) 
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TABLE VII. Values for the T mass in lattice units (sixth column) at the three different values 
of (3, derived from perturbative renormalisation parameters Z m and Eq (see text). 



P 




y 


aE 


aE NR 


a(M r ) ca i c 


5.7 


3.15 


1.25(6) 


0.45(20) 


0.5186(6) 


7.51(55) 


6.0 


1.71 


1.19(4) 


0.30(9) 


0.4537(5) 


3.91(23) 


6.2 


1.22 


1.31(10) 


0.21(4) 


0.3132(3) 


3.09(26) 
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